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Abstract 

Using thermodynamic and variational principles we examine a basic phase 
field model for a mixture of two incompressible fluids in strongly perforated 
domains. With the help of the multiple scale method with drift and our 
recently introduced splitting strategy for Ginzburg-Landau/Cahn-Hilliard- 
type equations [1], we rigorously derive an effective macroscopic phase field 
formulation under the assumption of periodic flow and a sufficiently large 
Peclet number. As for classical convection-diffusion problems, we obtain sys- 
tematically diffusion-dispersion relations (including Taylor- Aris-dispersion) . 
Our results also provide a convenient computational framework to macro- 
scopically track interfaces in porous media. In view of the well-known versa- 
tility of phase field models, our study proposes a promising model for many 
engineering and scientific applications such as multiphase fiows in porous 
media, microfluidics, and fuel cells. 

Keywords: diffusion-dispersion relations, homogenization, upscaling, 
multiscale expansion, porous media 



1. Introduction 

Fluid mixtures are ubiquitous in many scientific and engineering appli- 
cations. The dynamics of phase interfaces between fiuids plays a central 
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role in rheology and hydrodynamics [21 [3l HI El [6]. A recent attempt of 
a systematic extension towards non-equilibrium two-phase systems is j7], 
where the authors discuss the concept of local thermodynamic equilibrium 
of a Gibbs interface in order to relax the global thermodynamic equilibrium 
assumption. In [8], it is shown that the Cahn-Hilliard or diffuse interface for- 
mulation of a quasi-compressible binary fluid mixture allows for topological 
changes of the interface. 

The study of flows in porous media is a delicate multiscale problem. 
This is evident, for instance, by the fact, that the full problem without any 
approximations is not computationally tractable with the present computa- 
tional power. Also, from an empirical perspective the consideration of the 
full multiscale problem is very challenging due to the difficulty of obtaining 
detailed information about the pore geometry. These empirical and compu- 
tational restrictions strongly call for systematic and reliable approximations 
which capture the essential physics and elementary dynamic characteristics 
of the full problem in an averaged sense. A very common and intuitive 
strategy is volume averaging [9l [T0| . The method of moments [T2] and 
multiple scale expansions \13\ [TH [15] have been used in this context. The 
latter method is more systematic and reliable since it allows for a rigorous 
mathematical verification. The volume averaging strategy still lacks a con- 
sistent and generally accepted treatment of nonlinear terms. Therefore, the 
multiscale expansion strategy is used as a basis for the theoretical develop- 
ments in the present study. 

The celebrated works in |1H [T^ [TH] initiated an increasing interest in 
the understanding of hydrodynamic dispersion on the spreading of tracer 
particles transported by fiow, with numerous applications, from transport 
of contaminants in rivers to chromatogaphy. In [T7], it is shown that the 
multiscale expansion strategy allows to recover the dispersion relation found 
in [12] . The study of multiphase flows in porous media is considerably more 
complex; see e.g. the comprehensive review in [18] which still serves as a 
basis for several studies in the field. The central idea for many approaches 
to multiphase flows is to extend Darcy's law to multiple phases. With the 
help of Marie's averaging method [TH] and a diffuse interface model, effective 
two- phase flow equations are presented in \20\ [2T] . In [22] , Atkin and Craine 
even present constitutive theories for a binary mixture of fluids and a porous 
elastic solid. 

An application of increasing importance for a renewable energy infras- 
tructure are fuel cells [23j. This article combines the complex multiphase 
interactions with the help of the Cahn-Hilliard phase fleld method and a 
total free energy characterizing the fuel cell. An upscaling of the full ther- 
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modynamic model proposed in [23J is obviously very involved due to complex 
interactions over different scales. In this context, an upscaled macroscopic 
description of a simplified (i.e., no fluid flow and periodic catalyst layer) is 
derived in |24^ [25j . 

Consider the total energy density for an interface between two phases. 



e(x(X,t),t) := ^ 



dx(X, t) 



dt 



2 



^ I Vx0(x(X, t),t))f - ^F(<A(x(X, t),t)) , 

(1) 



where (p is a conserved order-parameter that evolves between different liquid 
phases represented as the minima of a homogeneous free energy F. The 
parameter A represents the surface tension effect, i.e. A oc (surface tension) x 
(capillary width). We will consider polynomial, homogeneous free energies 



F{<f>) := / f{s)ds, and /(s) := ass'^ + a2S^ + ais , 
Jo 



(2) 




The variable X stands for the Lagrangian (initial) material coordinate and 
x(X, t) represents the Eulerian (reference) coordinate. The last two terms 
in ([T]) form the well-known density of the Cahn-Hilliard/Ginzburg-Landau 
phase field formulation adapted to the flow map x(X, t) defined by 

^ - u(x(X,t),t) , 

The first term in ([T]) is the kinetic energy, which accounts for fluid flow of 
incompressible materials with the viscosity n, i.e., 

f t + (u-V)u-^Au + Vp = J7, 
\ divu = 0, 

where 77 is a driving force acting on the fluid. We are interested in the 
mixture of two incompressible and immiscible fluids of the same viscosity fi. 
Hence, we can employ generic free energies ([T]) showing a double-well form 
as is the case often in applications, e.g. [26j. 

Suppose that 17 C M'^, with d > the dimension of space, denotes the 
domain which is initially occupied by the fluid. Then, we can define for an 
arbitrary length of time T > the total energy by 

E{x):=[ [ e{x{X,t),t)dXdt. (5) 
Jo Jn 
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(Homogeneous case) < 



The energy ([s]) combines an action functional for the flow map x(X, t) and a 
free energy for the order parameter cj). This combination of mechanical and 
thermodynamic energies seems to go back to [271 EEl EHl El |8] . Subsequently, 
we will focus on quasi-stationary, i.e., Ut = and 77 7^ 0, and low- Reynolds 
number flows such that we can neglect the nonlinear term (u • V)u. Then, 
classical ideas from the calculus of variations [30j and the theory of gradi- 
ent flows together with an imposed boundary condition Jq^w{x) do{x) for 
u)(x) G H^/'^{d^) lead to the following set of equations 

—jiAu + Vp = rf in Qt , 

div u = in Qj- , 

u = on O^It ! 

g + Pe(u • V)0 = Adiv (V (/(0) - Ac/.)) in JIt , 
Vn(p '■= n • = w(x) on OQt , 

VnA(/) = on dQr , 

(/)(x, 0) = h{x) on O , 

(6) 

where Qt '■= J^x]0,T[, OCIt ■= i9^^^x]0,r[, A represents the elastic relax- 
ation time of the system, and the driving force r] accounts for the elastic 
energy [5] 

rf = -7div (V0 (g) V(p) , (7) 

where 7 corresponds to the surface tension [3T]. As in [32], we will set 
7 = A for simplicity. The dimensionless parameter Pe := ^"^^ is the Peclet 
number for a reference fluid velocity U := |u|, L is the characteristic length 
of the porous medium, and the diffusion constant D = krM obtained from 
the mobility via Einstein's relation for the temperature r and the Bolzmann 
constant k. We note that the immiscible flow equations can immediately be 
written for the full incompressible Navier-Stokes equations as in [5]. Our 
restriction to the Stokes equation is motivated here by the fact that such 
flows turn into Darcy's law in porous media |33| I34j . 

The main objective of our study is the derivation of effective macroscopic 
equations describing ^ in the case of perforated domains Q'' C M'^ instead 
of a homogeneous 17 C M'^. A useful and reasonable approach is to represent 
a porous medium = ^l'' L) B'' periodically with pore space 0,^ and solid 
phase B^. The arising interface dQ^ n dB'' is denoted by P. As usual, 
the dimensionless variable e > defines the heterogeneity e = j where 
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Reference cell Y 



Periodic covering by cells Y 



Homogenous approximation 




Figure 1: Left: Porous medium 57^ := Q\ as a periodic covering of reference cells 
Y := [0,1]''. Top: Definition of the reference cell F = U with £ = 1. Right: The 
"homogenization limit" e :— ^ scales the perforated domain such that perforations 
become invisible in the macroscale. 



£ represents the characteristic pore size and L is the characteristic length 
of the porous medium, see Figure [T] The porous medium is defined by a 
periodic coverage of a reference cell Y := [0, £i] x [0, ^2] x • • • x [0, (d], f-i G I^j 
i = 1, . . . , d, which represents a single, characteristic pore. The periodicity 
assumption allows, by passing to the limit e — )• (see Figure [T]) for the 
derivation effective macroscopic porous media equations. The pore and the 
solid phase of the medium are defined as usual by, 

Q.' := IJ e (y^ + z) n , := \J e [Y"^ + r\n = , 

where the subsets y^, Y^ C Y are defined such that VL'' is a connected 
set. More precisely, denotes the pore phase (e.g. liquid or gas phase in 
wetting problems) , see Figure [T] 

These definitions allow us to rewrite ^ by the following microscopic 
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formulation 



(Porous case) 



— e^/iAue + Vpe = - 


-7div (V(/)^ V(/)e) 


in fiy , 


div Ue = 




in Qfp , 


Ue = 




on , 


|</.e + Pe(u, • V)0, 


= Adiv (V - A<A.)) 


in nfp , 


V„(/), := n • V(/), = 




on /|, , 


V„A(/)e = 




on , 


M-,o) = H-) 




on rj^ , 

(9) 



where /|i := Px]0,T[. Even under the assumption of periodicity, the mi- 
croscopic system Q leads to a high-dimensional problem, since the space 
discretization parameter needs to be chosen to be much smaller than the 
characteristic size e of the heterogeneities of the porous structure, e.g. left- 
hand side of Figure [T| The homogenization method provides a systematic 
tool for reducing the intrinsic dimensional complexity by reliably averaging 
over the microscale represented by a single periodic reference pore Y. We 
note that the nonlinear problem ^ is characterized by a complex coupling 
between the micro- and the macroscale. This scale interaction leads to non- 
standard reference cell problems which also depend on the macroscale. As 
a consequence, the reference cell problems need to be computed for each 
macroscopic degree of freedom now. This high-dimensionality is an intrinsic 
feature of coupled non-linear partial differential equations studied in porous 
or highly heterogeneous structures, see [33 1211 1]- 

Obviously, the systematic and reliable derivation of practical, conve- 
nient, and low-dimensional approximations is the key to feasible numerics 
of problems posed in porous media and provides a basis for computationally 
efficient schemes. To this end, we relax the full microscopic formulation ^ 
further by restricting ^ to periodic fluid flow. By taking the stationary 
version of equation (|9|)4 on to a single periodic reference pore Y, we can 
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formulate the following periodic flow problem 



(Periodic flow) < 



— /iAyU + VyP = — 7diVy (Vyl/' (g) Vy^") IH , 

divy u = in , 

u = on dY"^ , 

u is y^-periodic , 
Pe(u • Vy)'0 = Adivy (Vy(/(V') - AyV')) in Y^ , 
Vnip '■= (n • Vy)'i/' = w on dY'^ , 

VnAyV- = on ay2 , 

tj; is y^-periodic . 

(10) 

We remark that in certain occasions it might be suitable to further reduce 
problem (10). For instance, if the reference cell is only fllled by one fluid 
phase, i.e., VyV' = 0, then one only needs to solve for the periodic Stokes 
problem (10)i-(10)4 under zero driving force, i.e., 77 = 0. It is also possible 
to impose a constant driving force 7/ := ei where ei denotes the canonical 
basis vector in the xi-direction of the Euclidean space. We further note that 



physically, the periodic fluid velocity defined by ( 10 ) can be considered as the 



spatially periodic velocity of a moving frame. Motivated by [MllIZ], we study 
the case of large Peclet number and consider the following distinguished 
limit. 

Assumption (LP): The Peclet number scales with respect to the char- 
acteristic pore size e > as follows: Pe ~ ^ . 

Let us first discuss Assumption (LP). If one introduces the microscopic 
Peclet number Pcmjc := ^77^) then it follows immediately that Pe = ^^^mic^ 
Since we introduced a periodic flow problem on the characteristic length 
scale £ > of the pores by problem (10), it is obvious that we have to apply 



the microscopic Peclet number in a corresponding microscopic formulation, 
see (11) below. Hence, the periodic fluid velocity u(x/e) := u(y) enters the 



microscopic phase field problem as follows 



(Microscopic problem) 



|'^e + ^(u(x/e).V),^, 
= Adiv(V (/(<Ae) - A</.,)) 
Vn(pe ■= n • V0e = w{x/e) 

VnAcj), = 

<^.(-,o) = v(-) 



in ilfp , 
on , 
on I|, , 
on 



(11) 



7 



We note that with our subsequently appHed upscahng strategy, we do not ac- 
count for boundary layers in the vicinity of rigid boundaries. Such boundary 
layers become increasingly important in the case of large Peclet numbers. f 
The main result of our study is the systematic derivation of upscaled im- 
miscible flow equations which effectively account for pore geometry starting 



from the microscopic system (lO)-(ll), i.e., 



(Upscaled equation) 



= div 
-C(x,0 



p/'(c/>o)M - [21^ - f{4>o)) - e(x, t) 
V0o) - /'(0o)div f (M, + K)V0o' 
*o)) , 

(12) 



+ ^div(M^V(div(DV(/)o 



where Q and C take the fluid convection into account. These two tensors 
account for the so-called diffusion-dispersion relations (e.g. Taylor-Aris- 



dispersion [HI US |T6] ) . The result (12) makes use of the recently proposed 
splitting strategy for the homogenization of fourth order problems in [Ij and 
an asymptotic multiscale expansion with drift (i.e., moving frame) intro- 
duced in [371 EH]. 

The manuscript is organized as follows. We present our main results in 
Section [2] and the corresponding proofs follow in the subsequent Section [3j 
Concluding remarks and open questions are offered in Section [4} 



2. Results: Effective immiscible flow equations in 

porous media 

For the presentation of our main result, we firt need to introduce the 
following: 

Assumption (LTE): (Local thermodynamic equilibrium) We assume 
that the following local thermodynamic equilibrium condition is satisfied in 
each reference cells of the periodic covering, i.e., 



5E{(P) 



^{(j)) = f{(f)) - x'^A(p = const. , 



(13) 



for each f = y part of the same cell Y , where fi stands for the chemical po- 
tential which is only allowed to vary over the different reference cells. 

The assumption (LTE) emerges as a key requirement for the homogeniza- 
tion of nonlinearly coupled partial differential equations. This equilibrium 
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requirement is well-known and generally accepted in physics if there exists 
a separation of macroscopic (size of the porous medium) and microscopic 
(characteristic pore size) length scales. This kind of equilibrium assump- 
tions are widely applied to a variety of physical situations such as diffusion 
[39], macroscale thermodynamics in porous media |40j . ionic transport in 
porous media based on dilute solution theory |35], and fuel cells [23], and 



density- functional theory |41t I42j. The local equilibrium assumption (13) 
emerges as a key requirements for the mathematical well-posedness of the 
corresponding cell problems which define effective transport coefficients in 
homogenized, nonlinear (and coupled) problems. 

As in U], the homogeneous free energy characterizing the equilibrium 
limiting values of two different phases needs to fulfill the following 

Assumption (F): Assume that the homogeneous free energy F satisfies 
for real parameters a2 > ai > 0, which define F as a double-well potential 
by F{s) = {s — ai)^(s — 02)^ , such that 

25(ai + - 20{al + al + 3aia2) > (ai + 02)^/4 . (14) 



This assumption is required for the well-posedness of the upscaled equa- 
tions. We summarize our main result in the following 

Theorem 2.1. (Effective convective Cahn-Hilliard equation) We assume 
that the Assumptions (LP), (LTE), and (F) hold and let G H'^{Q). 



Then, the microscopic equations (lO)-(ll) for immiscible flows in porous 



media admit the following effective macroscopic form after averaging over 
the microscale, i.e., 



pfiMM - - /'(0o)) - e(x, t) 

V(/>o) - /'('/'o)div f (M, + K)V^o 
Wo 

n • V0O = w{x/e) 



-C(x,t) 



+ ^div(M^V(div(DV</.o 



VnAcPo = 



where the tensors K := {Kki}Kk i<d 
D := {difc}i<i,fc<d, M^, = {m'^k}i<Lk<d' 



in 0,T 1 

on {do, \ I}rp , 
on {do, \ I}rp , 
in ft , 

(15) 



, e := {Okl}i<kJ<d' C := {cik} 



\<i,k<d' 



and Mu, = {m^^(x,t)} 



l<i,k<d 



are 
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defined by 



Kki 

Cik 

m" 



J_ x^d r _a_ 
\Y\ l^i,j=l Jyi dyi 



ac'''(y) 

\Y\ ^i,j=l Jyi dyi dvj 

P^/y,(u.Vy)C'='(y)dy, 



dy, 



ik 



\y\ 
1 



m. 



m 



^)*) |y| Z]j=i /yi f™jfc 



m, 



> dy, 
_ Pc,, 



(16) 



dy. 



The effective fluid velocity v is defined by v := -|^^ jyi'^'^ {v) dy where 
u solves the periodic reference cell problem The effective wetting 

boundary condition becomes :- 



Iy (aiXdYiiy) + a2XdY^{y)) dy. 
The corrector functions Q^^^^ G ^ir(^^) «^^^ ^ -^^^(f^; ^per(^^)) /' 



JL J_ 



1 < k,l < d solve in the distributional sense the following reference cell 
problems 



'■kl 



\-^d _d_ 
Z^i,j,k,l=l dyi ^ 

V„02 = 5e 

02 is Y -periodic . 



d /-kl 



on dY^ n dY^ , 



(17) 



Ed 
i,j,k=l 



dy, 
d 

k,i,j=l dyi 



Ed 
i.i. 



i,j,k=l ■ 



^ik 



V dy, 
\2 sr^d J)_ 
^ 2^k,i,j=l dy. 



mk 



/'(</.o)</'o™*J9%- 



/'(0o)0o™*J9%- 



,'^u;(y) ^■^ Y -periodic and Myii^w) — 0' 



Ed 

^d 



d 

=1 







Aviy) Y -periodic and AdyiiCv) 



in 



on 



on dY^ , 



(18) 



(19) 
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3. Proof of Theorem 



2.1 



As in [Tj, we introduce the differential operators 



^0 



E 



a 



Ed 
1,7 = 1 



d 
dxi 



+ 



E 



*J 8xj 

a 



ay, 



So 



E 



«.i=i ayi 



m, 



Ed 
1,7 = 1 



+ 



E 



dxi 
d 

V dxj 

a 

i,j=i axj 



m 



ij dyj 



(20) 



m. 



which make use of the micro-scale - =: y G y such that A^^ 



e + and B^ := e ^;Bo + e ^-^i + '^2. Herewith, the Laplace opera- 
tors A and div ^MV^ become Ati'^(x) = ^en(x, y) and div ^MV^ u"^(x) = 
i3e«(x, y), respectively, where n'^(x, t) := 'u(x — ^t, y,t). Due to the drift 
[571 EH] ) we additionally have 



d_ 



d_ 
di 



(21) 



where we find below by the Fredholm alternative (or a solvability constraint) 
that V := ^^FTj- Jyi ^Hy) dy . As in [1] we apply the method of formal asymp- 
totic multiscale expansions together with the splitting strategy introduced 
therein. As a consequence, we obtain the following sequence of problems by 
comparing terms of the same order in e, with the first three problems being. 



Oie- 



BqWo = Bofi4>o) - Pe„ 

no flux b.c. , 

wq is y^-periodic , 
AqVo = 

VnVo = 

(po is y^-periodic , 



(u-V)y)t;o inyi. 



in yi , 

on dY^ n ay^2 ^ 



(22) 
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Bowi = -Biwo - Bo [fiM'gl - Si/(</'o) 

-Pemic (u • Vy) vi - Pemic ((u - v) • V) vq in 

no flux b.c. , 

wi is y^-periodic , 
Aovi = -Awo 

^nVl = 

01 is y-'^-periodic , 
B0W2 = -A^ {B2W0 + Biwi) 



in , 

on dY^ n dY^ , 



(23) 



-Bo 
-Bi 



y"{cpo)4>'i+f'mH 



-^2/(<Ao) 

-Pemic (u • Vy) 02 - Pemic ((u - v) • V) Vi 

-dt{-A)-^wo 
no flux b.c. , 
W2 is y^-periodic , 
"4,0^2 = -A2V0 - Aivi + Wo 

V„f2 = Oe 

(p2 is y^-periodic , 



in y^ 



in yi , 

on ay 1 n dY^ , 



(24) 



As usual, the first problem (22) induces that the leading order terms 0o = 
+ 00 Wo are independent of the micro-scale y. The second problem 
(23) reads in explicit form for vi as follows, 



*,i=i dyi \S'->' 

= -div (efc - VyC^) = 

n • (Ve^ + efc) = 

,'^^'(y) is y-periodic and Alyi(C^) = 0, 



in yi , 

on dY}, n ay2 



(25) 



which represents the reference cell problem for vq after identifying vi = 
The cell problem for wi is substantially more involved since it depends 
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on the fluid velocity u and the the corrector from (25). Specifically, 

d 



E 



d 



... dyi\ dyj J dxk dyi 



m 



E 



m 



dxj 



E 



d 



d 



+ P-E f»o|-)«;g-Pe„.EK-v<)|| 



k,i,j=l 



i=l 



(26) 



where the last line can be simplified under the Assumption (LTE), i.e., 



Pe 

J 



2dw for 1 < A; < d, to 



(27) 



If we also rewrite the remaining terms under Assumption (LTE) as in [T], 
then we end up with the following cell problem, 

' _d_ 



I ^ ^mic \ 
"T f'C^^) Z^fc,i = l 



Pfi rn. J, n \ ^ I 



1 u' + v' 



f'{<f>o)4>o "-^dyj 



in 



J=l dy^ I 
+ /'(0o) 2^fe,i=l 1,9s/, U + V 

,C5;(y) is 1^-periodic and A^yi(^^) = 







on dY^ n ayj , 

(28) 



We are then left to study the last problem ( 24 ) arising by the asymptotic 



multiscale expansions. Problem (24)2 for V2 is classical and leads immedi- 
ately to the upscaled equation 



(29) 



ApUo := -div ( DVwo] = OiWq + go , 
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see also [T], where the porous media correction tensor D := {dik}i<,ik<d 
defined by 



Next, we apply the Fredholm alternative on equation (24)i, i.e. 



-Fe-mic (u • Vy) 02 



The term multiplied by in (31) can immediately be rewritten by 
' J - {B2W0 + Biwi) dy = -div (il^Vwo^ , 



where the effective tensor = is defined by 



m 



ik 



1 



The first two terms on the second line in (31) transform as in [T] to 



\Y\ 



-Bi 



B2fiM dy 



div 



9if'{4>o)M-2'-^M, 



where the tensor M^, = \ m^- > is defined by 

I ^ ) i<i,k<d 



ik 



dy. 



(30) 



{-A2 {B2W0 + Biw,) - Bo QricAo)-/-? + f'iM<p2 

fm^] - B2f{ct>o) - dt{-A)-'wo (31) 
Femic ((u - v) • V) fij dy = . 



(32) 



(33) 



(34) 



(35) 



The last term on the first line in (31) reveals the (non-trivial) coupling 
between the micro- and the macro-scale. After integrating over Y it can be 
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written in the following compact form, 
1 



1 



where 



-Bo 



:f"{Mri + f{ct>o^: 



dy = div (/'(</.o)M,V</)o 



/'((/>o)div (m,V</.o) +/'(0o)'^(x,t) 



(36) 



d_ 



m- 



(x,y,t) 



5%- 



(37) 



The terms in the last line of (31) become 

-TTF] I Pe^ic (u • Vy) 02 + Pernio ((u - v) • V) vi dy 
M I jyi 

= -TTTT / Pe^ic (u • Vy) (/)2 (iy 
1^ I jyi 



(38) 



Pe. 



1^1 



These considerations finally lead to the following upscaled phase field 
equation 



p— = div 



p/(0o)M - ( 2^ - /'(<Ao) ) - e(x,t) 



C(x,t) V0o) - /'(</'o)div ((M, + K)V</)o) (39) 



+— div ( MwV ( div ( DV^o ) - ^^o 



where the tensors K := {Kki}i<k,i<d^ © := {Gki}i<k,i<d^ and C := {cik}i<,^k<d 
are defined by 



^fcKx, t) := /^i (u • Vy)C'='(y) dy , 

C,fc(x,t) :=?^J^,(u'-y%,,^k^y)dy. 



dy, 



(40) 



The second order corrector C'^' is defined by (/)2(x, y, t) = Y2k i=i C^^(y) axt'^x;'' ' 
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where the variable (jP' is the solution of 



Vnh = ge on dY^ n dY^ , 

(p2 is -periodic , 



(41) 



which can be further simplified to 



v„02 = 5e on dY^ n 5yj 

i;^)2 is y^-periodic . 



(42) 



4. Conclusion 



The main new result here is the extension of the results in the study 
by Schmuck et al. in the absence of flow [1] to include a periodic fluid 
flow in the case of sufficiently large Peclet number. The resulting new ef- 



fective porous media approximation (15) of the microscopic Stokes-Cahn- 



Hilliard problem (lO)-(ll) reveals interesting physical characteristics such 



as diffusion-dispersion relations by (16)2^(16)3 for instance. The homoge 



nization methodology serves as a well-established tool for the systematic and 
rigorous derivation of effective macroscopic porous media equations starting 
with the fundamental works on Darcy's law [33} I34j . As a byproduct of the 
methodology developed here, one recovers rigorously and systematically the 
dispersion relations proposed in |12j . 

We note that the Cahn-Hilliard and related equations generally model 
more complex material transport \A2>\ |23] than classical Fickian diffusion. 
In this context, our result of an upscaled convective Cahn-Hilliard equation 
hence provides diffusion-dispersion relations for generalized non-Fickian ma- 
terial transport, that is, not just the product of a gradient of particle con- 
centration and a constant diffusion matrix. 

Of course, it would be of interest to remove the periodicity assumption on 



the fluid flow imposed by (10). This assumption implies a quasi-steady state 
on the fluid velocity and seems currently to be inevitable for the homoge- 
nization machinery to work such as the assumption of local thermodynamic 
equilibrium on the level of the reference cells. 
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So far, we are restricted to mixtures of two fluids by our model. An 
extension towards mixtures of > 2 components is studied in |44) where 
the authors compare a local and non-local model for incompressible fluids. 
In specific applications, it might be of interest to extend the here developed 
framework towards such multi-component mixtures. 

Finally, the rigorous and systematic derivation of effective macroscopic 
immiscible fiow equations provides a promising and convenient alternative in 
view of the broad applicability of the Cahn-Hilliard equations. The strength 
of our approach is based on its foundation on a thermodynamically moti- 
vated homogeneous free energy which is generally derived on systematic 
physical grounds. Moreover, even if a systematic derivation is impossible, 
one can often mathematically design such a free energy based on phys- 
ical principles and phenomenological observations. We expect promising 
insights of the systematically averaged interplay of micro- and macroscale 
in a wide variety of applications, e.g. for the optimization of microfluidic 
devices where the scale effects discovered here are crucial. A detailed nu- 
merical study might also allow to rationally analyze current applications in 
science and engineering and hopefully reveal potential new ones. 
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